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> Abstract 

2_j Recent LHC results suggest a standard model (SM)-like Higgs boson in the vicinity of 125 

If) GeV with no clear indications yet of physics beyond the SM. At the same time, the SM 

^^ is incomplete, since additional dynamics are required to accommodate cosmological dark 

(^ matter (DM). In this paper we show that interactions between weak scale DM and the Higgs 

^_J which are strong enough to yield a thermal relic abundance consistent with observation 

• • can easily destabilize the electroweak vacuum or drive the theory into a non-perturbative 

. !^ regime at a low scale. As a consequence, new physics — beyond the DM itself — must enter 

/\^ at a cutoff well below the Planck scale and in some cases as low as 0(10 - 1000 TeV) , a range 

jrt relevant to indirect probes of flavor and CP violation. In addition, this cutoff is correlated 

with the DM mass and scattering cross-section in a parameter space which will be probed 

experimentally in the near term. Specifically, we consider the SM plus additional spin 

or 1/2 states with singlet, triplet, or doublet electroweak quantum numbers and quartic or 

Yukawa couplings to the Higgs boson. We derive explicit expressions for the full two-loop 

RGEs and one-loop threshold corrections for these theories. 



1 Introduction 

With the quest for the elusive Higgs boson approaching a conclusive end, it is crucial that we 
evaluate the implications of its imminent discovery or exclusion for beyond the standard model 
(SM) physics. In the past, a conservative approach to this question has been to assess the 
consistency of the SM assuming a vast desert above the weak scale. Within this framework, 
numerous authors [1, 2, 3, 4] have analyzed the SM with respect to the stability of the elec- 
troweak symmetry breaking vacuum and the perturbativity of the underlying dynamics. As is 
well known, the relevant physics is determined by running couplings into the ultraviolet using 
renormalization group equations (RGEs). Hence, results depend sensitively on the weak scale 
boundary conditions, among which the mass of the Higgs boson is perhaps most critical. Recent 
experimental results from ATLAS [5] and CMS [6] suggest a value of the Higgs mass around 125 
- 126 GeV. This value corresponds to the SM with perturbative couplings up to the Planck scale, 
and a metastable electroweak vacuum with a lifetime longer than the age of the Universe [2, 3]. 

Despite its successes, however, the SM is almost certainly incomplete, since new physics is 
required below the Planck scale in order to accommodate observed phenomena such as dark 
matter, baryogenesis, neutrino masses, and the QCD theta parameter. In principle, the new 
dynamics may couple directly to the Higgs boson, therefore inducing deviations from the usual 
vacuum stability and perturbativity bounds of the SM. At the same time, many of these theories, 
e.g. the neutrino seesaw and the QCD axion, are accommodated by rather high scale dynamics, 
in which case there will be a minimal effect on the running of couplings. A notable exception 
to this is weakly interacting massive particle DM, whose mass scale is constrained to be low by 
thermal freeze-out. 

In this paper we carry out a general analysis of vacuum stability and perturbativity in the SM 
augmented by weak scale DM. We consider additional spin or 1/2 states with singlet, triplet, 
or doublet electroweak quantum numbers, and including all quartic and Yukawa couplings to 
the Higgs boson consistent with the stability of DM. For our analysis we have derived two-loop 
RGEs and one-loop threshold corrections for these theories. We find that fermionic DM tends to 
destabilize the electroweak vacuum, such that new dynamics below the Planck scale is required 
ensure stability. On the other hand, scalar DM tends to stabilize the vacuum, though demanding 
perturbativity by itself may require new physics at intermediate scales. Moreover, imposing a 
thermal relic abundance fixes a minimum value for the new quartic or Yukawa couplings which 
can have a substantial effect on vacuum stability and perturbativity. This effect is correlated 
with the DM mass and scattering cross-section in a parameter space which will be accessed by 
DM direct detection experiments in the near future. In many cases the required cutoffs for these 
theories can be as low as 10 - 1000 TeV, a range of energy scales which is reachable with indirect 
searches, such as those probing flavor and CP violating interactions. Our results are consistent 
with previous analyses of specific scalar [7, 8] and fermionic [9] DM models. 

Let us elaborate briefly on the notion of vacuum stability in the SM and beyond. Depending 
on the Higgs quartic Xh at the scale A, our vacuum may be any of the following: 

i) Stable {Xh > 0). The vacuum is the absolute minimum and will never decay. 
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Figure 1: Vacuum structure of the SM as a function of the Higgs boson mass. Regions of 
stabihty/metastabihty/instabihty are denoted in blue/purple/red respectively. The solid lines 
indicate central values while the dotted lines indicate ±2o" error bars on the experimental mea- 
surement of the top quark mass. 



ii) Metastable (0 > Xh > A^). The vacuum is not the absolute minimum, but its lifetime is 
longer than the age of the Universe. 

iii) Unstable (Xh > Xh)- The vacuum is not the absolute minimum, and it decays within the 
age of the Universe. 

Here the critical coupling Xh is determined by the requirement that the tunneling rate per unit 
volume is comparable to the age of the Universe. In particular, we demand that if ^ = F, where 
H^^ c^ 3.7 Gyr and F reads. 



max 



R-^exp{-167T^/3\XH\) 



R-^<A 



Here R is the characteristic length scale of the bounce, which is bounded by the cutoff. As 
we will elaborate on later, the vacuum structure may be more complicated if the new physics 
includes additional scalar particles. 

It is possible that our vacuum resides in a stable or metastable regime, but the unstable 
regime is of course excluded by our existence. In much of our analyses, it will be convenient to 
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Figure 2: Metastability bands for SM + singlet /doublet and triplet /doublet fermion, shown as 
a function of the Higgs mass. Regions above/below each band are stable/unstable. Each band 
is labeled with the corresponding value of the Yukawa coupling, i/s^t = Vst- The dashed lines 
correspond to the value of ttih suggested by [5, 6]. 



summarize the nature of the vacuum by depicting the "metastability band" in parameter space 
defined by the region 



Xh < Xh < 0. 



(2) 



Note that for theories in which supersymmetric dynamics enters at the cutoff A, i.e. split or 
high-scale supersymmetry, absolute stability is required by the fact that the potential is positive 
semi-definite in all field directions, so Xh > 0. 

By running the RGEs into the ultraviolet, subject to infrared boundary conditions, it is 
possible to determine the scale A at which the quartic coupling crosses through the metasta- 
bility band. For example, in Fig. (1) we plot the scale A indicating the onset of stabil- 
ity /metastability /instabilility in the SM as a function of the tuh- Our results are in nice 
agreement with existing calculations in the literature [2, 3, 4]. 

The remainder of this paper is as follows. In Sec. 2 we briefly summarize the set of models 
to be studied and establish notational conventions. We present our methodology in Sec. 3 
regarding the running of couplings and evaluation of DM properties. Finally, in Sec. 4 we 
discuss our results. The two-loop RGEs and one-loop threshold corrections for these theories 
are presented in App. A and App. B. 
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Figure 3: Higgs mass bounds as a function of the cutoff, for the SM + singlet scalar. The left 
panel depicts metastability bands and the right panel depicts the scale at which the largest 
coupling in the theory becomes non-perturbative. Each label denotes the corresponding value 
of the cross-quartic coupling, ks, and we have fixed the self-quartic coupling, A^ = 0. 



2 Models 

In this section we briefly summarize the models we will analyze and provide notational conven- 
tions. We choose the following normalization for the SM quartic and Yukawa couplings. 



C = ytqeHe + y,qH^h- + yJH^T' + ^{\H\^-v^)\ 



(3) 



where v = 174.1 GeV is the Higgs vacuum expectation value. Throughout, we ignore the effects 
of the light fermion generations because they are negligible for our results. 

We will augment the SM with new states, sending £ — t- £ + AC. Throughout, if a particle 
X is a fermion, then yx will denote its Yukawa coupling of X to the Higgs. If a particle X 
is a scalar, then Ax will denote its self-quartic coupling, while kx will denote its cross-quartic 
coupling to the Higgs. Finally, the mass parameter for the particle X will be denoted by mx- 
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Figure 4: Same as Fig. (3) but for SM + triplet scalar and varying kt with At=0. 



For the case of new ferniions we consider the following theories: 



singlet/doublet fermion: - AC = husS'^ + niDDD' + ysHSD + y'gH'SD'' 
triplet/doublet fermion: - A£ = -mxT"^ + moDD^ + yrHTD + y^H^TD'', 



(4) 



where we have defined H'^ = eH*, and we have included all renormalizable operators permitted 
by gauge symmetry. Here S and T are Majorana fermions and D and D'^ form a Dirac fermion. 
As a consequence of the additional Yukawa couplings to the Higgs, the new fermionic states will 
generally mix after electroweak symmetry breaking, although there is a preserved Z2 symmetry 
which makes the lightest among these fields a DM candidate. We do not consider the fully 
mixed singlet /triplet /doublet theory, since all the qualitative features of this model are already 
evident in the singlet /doublet and triplet/doublet cases, while the proliferation of parameters 
would make results difficult to present. 

The fermionic theories above are of course a generalization of the bino, wino, and Higgsino 
sector of the MSSM. In the chosen parameterization, the precise correspondence between these 
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Figure 5: Same as Fig. (3) but for SM + doublet scalar and varying ko with Ad = A'^ = fi:'^=0. 



theories and the MSSM is 

S ^ B 
T ^ eWe^ 
D ^ eHd 



ms i^ Ml 
rriT ^ M2 



ys O 9'JV2 

Vt ^ -9d/V2 
Vt ^ 9u/\/2 



using the notation of [10]. Furthermore, in the limit of exact supersymmetry, gi 
and gj^ = g'^'^ cos j3. 

For the case of new scalar s we consider the following theories: 
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singlet scalar: — AC 
triplet scalar: — AC 
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doublet scalar: - AC = niDW + ^\Dr + ^IDHH 



DH^?. 



(6) 



Here S and T are real scalars while D is a complex scalar. Contrary to the fermionic case, the 
pure doublet scalar case can have direct couplings to the Higgs and therefore can be considered 



alone without including singlets or triplets. For the singlet and triplet theories we have included 
all operators permitted by gauge symmetry subject to a discrete Z,2 symmetry under which 5* 
and T are odd. For the doublet scalar theory we have assumed a Peccei-Quinn symmetry which 
acts oppositely on D and H in order to reduce the number of possible operators to a manageable 
number. Hence, all of these theories carry a Z2 symmetry under which the lightest odd particle 
is a prospective DM particle. 

3 Methodology 

To study vacuum stability and perturbativity in the presence of additional dynamics we have 
produced code which inputs an arbitrary Lagrangian and outputs the corresponding two-loop 
RGBs and one-loop weak scale threshold corrections to Higgs and top couplings. For the two-loop 
RGEs we employed the MS expressions of [11], including the corrections discussed in [12]. We 
have also computed one-loop weak scale threshold effects for the theories under consideration 
using the methodology of [13]. All of these results, including the general formulae for the 
threshold corrections using the notation of [12] and [14], are presented in detail in App. A and 
App. B. 

As is well-known, the running of \h is highly sensitive to the weak scale values of A^, yt and 
as- This is because Higgs boson/top quark loops are the primary contribution driving A^ to be 
positive/negative in the ultraviolet. In turn, the strong interactions feed into the top Yukawa 
coupling. In our analysis we take the current values of the top quark pole mass and the MS 
running QCD coupling measured at the Z pole [15]: 

mt = 173.2 ± 0.9 GeV (la) (7) 

asimz) = 0.1184 ±0.0007 (la). 

On the other hand the Higgs boson has not been found and its mass has yet to be measured. 
Nevertheless, we make use of recent experimental results from the LHC. Since almost all of the 
theoretically allowed region for ttlh is excluded at 95% CL by both ATLAS and CMS, and both 
collaborations observe an excess at 125 - 126 GeV, assuming that the present excess corresponds 
to the Higgs boson, we will use 

niH = 125 GeV, (8) 

in our analysis. We relate the top quark and Higgs boson pole masses to yt{Q) and Xh{Q), the 
Yukawa and quartic couplings renormalized at an MS scale Q, via 

mt = yt{Q)v[l + 6t{Q)]~' (9) 

mjj = 2XiQy[l + 6Hm- (10) 

Here 6t^H includes one-loop threshold corrections from the SM and new physics contributions. 
We present general formulae for the one-loop threshold corrections in App. B. Furthermore 6t 



includes the shift between pole and running mass due to QCD interactions which is known to 
three loops [17]. 

We have chosen our MS matching scale Q to be the top pole mass rut. Our boundary 
conditions for SM parameters are determined by Eq. (8), Eq. (8), and Eq. (10), where we have 
run as{rnz) from niz to rrit using the SM two-loop RGEs enhanced by the three-loop QCD beta 
function [15]. We then run the RGEs up to a scale Q = A and determine the value of Xh at 
that scale. As defined in Eq. (1), ii X^ < Xh < ^ then the vacuum is metastable but long-lived, 
while if Xh < Xh then the vacuum is unstable and decays within the age of the Universe. 

Note that Xh > is a necessary but not sufficient condition for absolute stability of the 
vacuum. Likewise, > Xh > Xh does not guarantee that the hfetime of the vacuum exceeds the 
age of the Universe. The reason for this is that in theories with new scalar fields, the vacuum 
structure is enriched. Absolute stability of the vacuum implies further conditions on scalar field 
theories: the self-quartic couplings satisfy A^, At, Xd > 0, while the cross-quartic couplings are 
bounded by 



f^s > -2^XhXs (11) 

KT > -2^/XhXt (12) 

^d + ^'d > -I^XhXd, (13) 

for all scales below the cutoff A. If these criteria are not satisfied, then there will exist field 
directions in which the potential is unbounded from below at large field values. Alternatively, if 
the additional scalar fields acquire vacuum expectation values, this can also substantially alter 
the stability of the vacuum [16]. However, we will not consider this possibility since our interest 
is in DM. 

In addition to the question of stability we also investigate the perturbativity of interactions 
in the ultraviolet. For our purposes we define perturbativity to be the criterion that for each 
coupling (yf, the contribution of g to its own beta function is bounded by unity. Precisely, we 
require that dg/dlogQ = (5{g) < 1, which can be read off trivially from RGEs presented in 
App. A. As we will see, the constraint of perturbativity will be largely unimportant for the case 
of new fermionic states, but can play an essential role in determining the cutoff for theories with 
new scalars. For the scalars, the perturbativity bounds on the scalar couplings are 

4 < 167rV4, A| < 167rV36 (14) 

4 < 167rV4, A| < 167rV44 (15) 

4,4^<167rV2, A^ < 167rVl2. (16) 

To evaluate the properties of DM in each of these theories, we have implemented each 
model in LanHEP [18] and evaluated relic abundances and direct detection cross-sections in 
micrOMEGAs [19]. For the direct detection cross-sections we have employed the most recent 
lattice results [20] for the nuclear form factors, 

f!^l = 0.0280 4^j = 0.0280 f^^] = 0.0689. (17) 



Nuclear uncertainties will in general affect these results, particularly for the strange quark 
contribution. For our direct detection bounds, we have taken the current limits derived from 
the XENONIOO experiment [21]. 

Lastly, let us briefly comment on electroweak constraints. For the range of allowed couplings, 
the theories described in Sec. 2 typically possess violations of custodial symmetry and induce 
corrections to the p parameter. However, for all of the theories we consider here, we have checked 
that these effects are well within experimental bounds. We also have checked the corrections to 
the S parameters across the whole range of parameter space considered here and found them to 
be small as well, mostly due to the relatively heavy spectra. 

4 Results 

In this section we present the results of our analysis. We will begin with a general discussion 
of stability bounds, independent of the fact that the new states may constitute the DM of 
the Universe. Afterwards, we will present an analysis dedicated to DM in which we impose 
considerations from thermal relic abundance and direct detection constraints. 

4.1 Stability, Perturbativity, and New Physics 

We now analyze the effect of new physics on the vacuum as a function of the model parameters 
input at the weak scale. Consider Fig. (2), which characterizes the stability of the vacuum 
in the (A, 171^) plane for the singlet /doublet and triplet/doublet theories outlined in Sec. 2. 
These plots are similar to the SM plot presented in Fig. (1). Each colored region corresponds 
to the metastability band, A/^ < A^ < 0, for a different value of the Yukawa coupling, fixing for 
simplicity ys = y% and yx = yx- Hence, the upper and lower boundary of each band denotes 
the scale A at which A// = and \h = Xh, respectively. In this plot we have fixed ms = 100 
GeV, rriT = 200 GeV, and niD = 500 GeV, entering the one-loop threshold corrections. 

Fig. (2) demonstrates how new Yukawa couplings to the Higgs tend to destabilize the vac- 
uum. Indeed, the scale of instability drops precipitously — to within even a couple of orders of 
magnitude of the weak scale for uih — 125 GeV — as soon as y^lp > 0.5. The reason for this 
trend is that loops of new fermions tend to drive the Higgs quartic negative in the ultraviolet. 
To give a sense of the relative size of these couplings, recall that in the MSSM at tan/3 = 1, 
the analogous gaugino couplings correspond to ys = Vs = 5'V2 — 0.2 and yr = yr = 5'/2 — 0.3. 
Lastly, note that the constraint of perturbativity is typically unimportant. In particular, current 
Higgs boson exclusions [5, 6] suggest that rriH is in the light range, so \h is small and remains 
perturbative. Likewise, as long as y^j^ are not too large at the weak scale, couplings do not run 
non-pert urbative. 

In the left panels of Fig. (3), Fig. (4), and Fig. (5) we present metastability bands in the 
(A,m//) plane for the singlet, triplet, and doublet scalar theories, respectively. Here we have 
fixed ms = 100 GeV, tut = 200 GeV, and rriD = 500 GeV for the one-loop threshold corrections. 
Meanwhile, we have assumed vanishing self-quartic couplings at the weak scale, A^ = A^- = A^ = 
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Figure 6: Contour plots for SM -|- singlet/doublet fermion with singlet-like DM, shown in the 
{yg, ys) plane at fixed values of ms and mo- The purple bands corresponds to flh'^ = 0.11 ±0.01. 
The red/blue regions are excluded/allowed by XENONIOO, with asi denoted. The gray contours 
in the upper left/lower right quadrants denote the scale A (GeV) at which the vacuum becomes 
metastable/unstable. 
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Figure 7: Same as Fig. (6) but for SM + singlet /doublet fermion with doublet-like DM. 
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Figure 8: Same as Fig. (6) but for SM -|- triplet /doublet fermion with triplet-like DM. 
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Figure 9: Same as Fig. (6) but for SM + triplet /doublet fermion with doublet-like DM. 
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A'^, = 0, whose effects on Higgs stability are subleading, while the cross-quartic couplings have 
been allowed to vary. For the doublet scalar theory we have also fixed k'^, = for simplicity. 
Regions corresponding to an unstable vacuum, i.e. Xh < Xh, are not shown because these occur 
in a light Higgs range which is already experimentally excluded. 

From these plots we see that additional cross-quartic couplings to the Higgs tend to stabilize 
the vacuum. This occurs because scalar loops drive the Higgs quartic positive in the ultraviolet. 
Thus, while Fig. (1) indicates that tuh — 125 GeV implies a stability cutoff of A ~ 10^° GeV, this 
can be lifted up to the GUT or Planck scale, A ~ 10^^ — 10^^ GeV, if the cross-quartic couplings 
are evenly modestly sized. Finally, note that in these plots all of the theory parameters remain 
fully perturbative, primarily because the weak scale values of the cross-quartic couphngs were 
chosen to be small. 

On the other hand, the right panels of Fig. (3), Fig. (4), and Fig. (5) depict the scale A at 
which any one of the theory parameters become non-perturbative. We use the simple criterion 
for perturbativity outhned in Sec. 3. We see that perturbativity by itself places a significant 
constraint on the scale where new physics must enter, especially for larger couplings. 

4.2 Stability, Perturbativity and Dark Matter 

Thus far we have neglected the fact that the states under consideration should comprise a viable 
thermal relic DM candidate consistent with current direct detection constraints. Imposing these 
additional criteria reduces the theory parameter space substantially, allowing for more conclusive 
statements. 

To begin, let us consider the case of the singlet/doublet fermion theory represented in Fig. (6) 
in the {y%,ys) plane. Different values for irts and mr, are shown in different plots. Because we 
have chosen ms < ttlj:,, the DM is singlet-like. Each point in this plot corresponds to a different 
choice of model parameters, for which we can evaluate stability and check consistency with 
respect to a thermal DM relic abundance and current constraints on asi, the spin independent 
DM scattering cross-section, from XENONIOO. The red/blue regions are excluded/allowed by 
XENONIOO, while the purple bands correspond to the WMAP favored region, Qh'^ = 0.11 ±0.01. 
As denoted by the key, the gradations of blue denote different values for the spin independent 
direct detection cross-section, asi- The labeled contours in the upper left/lower right triangles 
denote the lower/upper values of A which border the metastability band. Said another way, 
the contours in the upper left triangle correspond to A^ = while those in the lower right 
correspond to Xh = Xh- 

The shape of the XENONIOO excluded regions in red can be understood as follows. As 
the Yukawas increase, the coupling between the Higgs boson and DM tends to increase as 
well, modulo the effects of mixing. Consequently the regions excluded from constraints on asi 
typically occur at larger values oii/g . An important caveat, however, is that the Higgs coupling 
to DM is strongly suppressed for certain values of yslv'si since they control the mixing angles. 
This effect is responsible of the allowed regions in blue extending to the boundaries of the plot. 
This feature was discussed in the case of general mixed singlet /doublet DM in [22]. 

Likewise, the shape of the WMAP favored region is clear because at ys = '^y's-, the coupling 
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Figure 10: Contour plots for SM + singlet /doublet fermion with singlet-like DM, shown in the 
(mo,m+ — mo) plane at fixed values of y%/ys- All points are consistent with Vth"^ = 0.11, while 
the red/blue regions are excluded/allowed by XENONIOO, with asi denoted. The solid gray 
contours denote the scale A (GeV) at which the vacuum becomes metastable. 
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Figure 11: Plots for SM + singlet scalar DM, shown as a function of the quartic couplings, 
Ks and Xs- The left panel depicts contours A (GeV) required by stability and perturbativity. 
The red/blue/purple regions correspond to regimes in which A is determined by metastability 
due to a negative Higgs quartic/metastability due to a large and negative cross-quartic/non- 
perturbative couplings. The right panel depicts the value of rris required for a thermal relic 
abundance along with the predicted asi- 
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Figure 12: Same as Fig. (11) but for SM + triplet scalar DM. 
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Figure 13: Same as Fig. (11) but for SM + doublet scalar DM. 



between the Z boson and DM vanishes. Hence, DM annihilation channels involving the Z boson 
are closed. In order to boost the annihilation cross-section, the Yukawa couplings must then 
increase, enhancing annihilations to Higgs bosons. 

Finally, let us address the contours characterizing vacuum stability. As expected, as ys and 
t/g increase, the Higgs quartic is driven negative in the ultraviolet and so the vacuum grows more 
and more unstable. Because metastability is a less stringent criterion than absolute stability, the 
curves in the lower right triangle are in general larger than those in the upper left triangle. From 
the plots it is clear that fermionic DM produced from thermal freeze-out and consistent with 
XENONIOO tends to destabilize the electroweak symmetry breaking vacuum. In particular, the 
vacuum becomes metastable in some cases even below A < 10^ GeV. The fact that the cutoff A 
of this theory may be so low indicates that in these models new physics may be lurking in an 
energy range relevant for indirect searches in flavor or CP violating interactions. Note that if 
the stable fermion is only a subdominant component of the DM in the Universe, then the purple 
band is a lower bound on the couplings ys and y^, in which case the vacuum will be destabilized 
even more. 

Fig. (7) is the same as Fig. (6) except mo < ms, so DM is doublet-like. Likewise, Fig. (8) 
and Fig. (9) are the same as Fig. (6) only for the doublet /triplet fermion theory. In Fig. (8), 
we have chosen rriT < ttlj:, to yield triplet-like DM, while Fig. (9) corresponds to m£, < rriT and 
thus doublet-like DM. 

Next, consider Fig. (10), in which we present contour plots in the (mo, "01+ — itiq) plane for 
the singlet /doublet fermion theory. Each panel corresponds to a different value of the mixing 
angle, y^/ys, and we have again fixed ms < mo so that DM is singlet-like. Every point on 



this plot is consistent with the WMAP favored region, Qh"^ = 0.11, thus fixing one combination 
of the four theory parameters {ms.mo.ys^y's}- We have then plotted contours of the scale A 
at which the vacuum becomes metastable, as well as contours of a si- As before, the red/blue 
regions indicate parameter space which is excluded/allowed by XENONIOO. 

Some remarks are in order. As we saw in Fig. (6), at certain values of y%/ysi there is 
a cancellation in the coupling of the Higgs boson to DM due to mixing, and asi drops to 
zero. In practice there will still be higher order contributions to asi-, but they will be highly 
suppressed and therefore we neglected to present them here. On the other hand, in these regions 
of parameter space the spin-dependent cross-section, asD will be non-zero, since it is controlled 
by Z boson exchange. In particular, aso lies in a range of 10~^^ - 10~^^ cm^, accessible to 
experiments in the upcoming future. The same feature appears in Fig. (10), where the white 
regions corresponds to the point of maximal cancellation in the direct detection scattering cross- 
section. From this plot, it is clear that the next generation of direct detection experiments will 
limit a large portion of the allowed region. 

Regarding the stability bounds, one sees that as r7i+ — uiq increases in value, then the singlet 
and doublet fermion components in general become less mixed. In order that the singlet-like 
DM not be overly abundant, ys and y^ must increase, which in turn destabilizes the vacuum. 
As can be seen from Fig. (10), the vacuum can become metastable at scales as low as 10^ GeV 
in the case of large ttlj^ — itlq. 

Lastly, let us now consider the case of scalar DM. We begin with Fig. (11), in which we 
present results for the SM + singlet scalar DM. In the left panel of Fig. (11), each colored 
region depicts a distinct phase of the theory. For any given point in parameter space, a cutoff A 
may be required due to some pathology — either because a stability condition fails for the Higgs 
quartic {\h < 0), a stability condition fails for the cross-quartic {ks < —2\/^XhXs), or a coupling 
becomes non-perturbative. These regimes are indicated in the left panel by the blue, red, and 
purple regions, respectively. Meanwhile, the labeled contours indicate the lowest cutoff dictated 
by any of these criteria failing. Naturally, the stability bound on the cross-quartic dominates for 
negative values of ks, while the perturbativity bound dominates for large values of A5. In the 
right panel of Fig. (11) we have plotted the critical value for ms required to yield the WMAP 
favored relic abundance, Qh"^ = 0.11, as function of the cross-quartic, ks- These values for ms 
were used to compute the one-loop threshold corrections in Fig. (11). Fig. (12) and Fig. (13) 
are the same, but for SM + triplet and doublet DM. 

5 Conclusions 

In this work we have studied the effect of weak scale DM on vacuum stability and perturbativity. 
Our primary result is that theories with thermal relic DM coupling via the Higgs typically require 
new dynamics in addition to DM at a scale A well below the Planck scale. For heavy fermionic 
DM, A can easily be as low as 10 - 1000 TeV as a consequence of vacuum stability. On the 
other hand, in theories of scalar DM, A takes a broad range, depending largely on the criteria 
of perturbativity as well as stability. Because A is correlated with the mass and scattering 
cross-section of DM, direct detection experiments may serve as an indirect probe of the scale 
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of new physics. In short, while new physics at the scale A is likely inaccessible to the LHC, 
indirect effects may be probed in the cosmology frontier via direct detection experiments, or in 
the intensity frontier via high luminosity probes. Lastly, note that the new dynamics introduced 
is often required to couple to the Higgs — e.g. if its purpose is to prevent vacuum decay — but it 
may or may not interact directly with SM fermions. In such cases, experiments probing flavor 
conserving CP violation may be more likely venues for indirect signals than those considering 
flavor violating phenomena. 
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A Two-Loop Renormalization Group Equations 

We define the following convenient notation in which the RGEs for a given coupling strength g 
are given by 

da 4" , of 



+ 



dt (47r)2 (47r)4' ' ' 

where t = log (Q/Qo) and bg and bg are the one- and two-loop beta functions for g, respectively. 
For thoroughness we also recall the well-known RGEs for the SM, defined in the notation of 
Eq. (3). Afterwards, we present results for new physics models, which induce contributions of 
the form 

^9^ ^ &i'^ + A6(i) (19) 

fef ^ 6« + A6W, (20) 

for each coupling g. When the expression Abg is not shown, then the corresponding one- or 
two-loop quantity receives no corrections from new physics. The RGEs for the fermionic models 
were derived independently and presented in [4]. 
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A.l SM 



3 



, (1) _ 41gf 
"91 — 10 

/,(2) _ _1„3 2 _ 17 3 2 _ 3 3 2 , 199^1 , 27 2 3 , 44„2 3 



(1) _ 19sf 



^92 



4? = -Uvl - Uv't - blvl + ¥ + Mai + 12^1^1 
by^ = Ivlvt - Egbt ~ hht - ^glvt + ytVr + — 



WbVt - ^gtvt - Igm - mvt + ytVr + -f- 

by} = iogbbVt + Mvlvt + "^glylvt + Ivlytvl - ^vlvl - Ivtyt + ^glvtvl + ^glvtvl + ^glv^ + 
f ^ii/*^ + ^QgM + H^f y* - f ^ly. - i^Hvt - Maht + Mglvt + 9^|^ii/, - eA^y? + 



byj = -\glyb - Iglvb - ^glvb + Ivbv'i + Vbvl + ^^' 



4? = mglvby^t + f^glvbyl + W^VbV^t + f ^?i/fei/' + 'iglvbvl + ff ^^y^' + ''^glvl + se^^fy^ _ 121^4^^ 
fg% - iQSgivb - ^ghhb + %glglyb + '^glglvb - QylXn + Ivb^l + Ivbvhl - ^ivlvl 

Ivhvf — T-u?-u? — ^Vhvt — 12-?;? 



Ivbyf - Ivbr - Ivbvt - 12l/, 
€ = Hvr - Iglvr - Iglvr + Sy^. + ^ 



ighlyr+fglyhr+^Oglybr+lybtyr-^ybyr - ^ytyr+'igbtyr+fghtyr+'^^gl 

537 ^2„,3_|_ 165 ^2„,3_i_ 1371 ^4„, 23_4„, i 27_2„2„, c\ „,3 i 3 ^2 „, 27„,2„,3 27„,4„, 



6i^] = 121/2 A^ - 121/4 - IgfX^ - ^glXu + If + Mai + ¥ + 12AHt/,2 + 4Ah1/2 + 12A^ - I2yt - 4^/^ 



i^?y,^ Ah + 'iglylXn + 80(/|i/2 Ah + l^g\yl + 1^2^,4 + f ^2^2^2 _ 64^2^4 _ |^4^2 _ 42yl\Hyl - 
72yl\l - ?,yt\H - I2ylyt - I2yty^ + QOyt + 'igjXny^ + fglXny^ + SOglXny^ + fgfXnyl + 
fgl^Hyl+'^gfXH+fglXl+^^glgjXH+^AglxWigl^H-'i^gtyrf^^^^^ 

RAfj^oA - 9„4 2 _ 9 4 2 _ 24 2 4 , 33 2 2 2 _ 3 4 2 _ 3411gf _ 1677 2 4 _ 289 4 2 , 305g| _ 

^^gd,yt 2"2yt 2"ii'T 5i/iyrT^ 5 y2yiyT 2y2yr 1000 200 "2yi 40 y2yi T^ 8 
72A2,i/2 _ 3A^^4 _ 24A2^^2 _ ^^^4 _ 73^3^ + gg^e + 2O1/6 
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A.2 SM + Real Singlet Scalar 

A4f = \4yt 



A/,(2) _ 12 

AC = 4 



A6f] = -5A^k| - 44 

C = 4 + 36A| 

C = -61/64 + l^f 4 + %i4 - 64l/' - 241/? - 204As - 44 - 816AI 

&S = %b/«5 - jo9ff^s - i^i/«5 + 6AH/t5 + Qf^sVt + '^f^sVr + ^'2ks\s + 44 

SQglXnus + ^gl^sy^ + fgi^sV^ + ^Ogl^sy^ + ^gf^syl + fg^Ksy^ + 'S^gtKs + IgWs + 
Ighli^s + 'igWs - ^-§gt^^s - SQXn^^sy^ - l2\HKsyl - IhXlKs - 36\h4 - 124t/2 - 

f^syt - ^4y'r - l^sy'r - 120/.5AI - 724A5 - ^ 



A.3 SM + Real Triplet Scalar. 



A ^(1) _ 9l 



^c 



A6£ 

a4' 

A6i^ 



A , (2) 289S 

= gtyt + |i/t4 
= gtyh + |i/b4 

— 6*2 l/r + 4^f^TyT 

A6g = f (^lA^ + 3(}glKT + ASgWr - ¥ " n)^?^2 " 15Ai,4 - 124 
6i'j = -2AglXT + Qgl + 4 + 44A| 

C = -%^T + 10^2V + 6^i4 + |^?4 + f ^|AT + 512^iA|-ii^-6t/,24-24y^-204AT- 
44 - 1104A^ 

C = ^ghW + l^fy^T + ^QglyW - Qgtyl - sey^A^/t^ - 21^2^2^^ - 12^,^4 - Wb^T + 
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2>Qgl\HHT + ^igl\H^^T + SO^^^A^ + ^gly^KT + ^igWtKT + mgly^^Kr - Gg^y^ + fglnryl + 
fgJKry', + 320glKTXT + '-§glt^T + HgWr + IqIqW + |^?4 + ^^dW + IOO^^^'At - 
'igbl + ^ - Iglgt - SQXny^KT - VlXuK^t - ^^X\Kr - 36A^4 - 11ylK\ - fyf^r - 
Anlyl - l^ryt - 200/tTA^ - 1204At - ^ 



A. 4 SM + Complex Doublet Scalar 



Ab^' 



'91 



Al^' 



'91 



Abg 

A6g 
Afcg 



A6 



.3 



= ^ 

10 



-Ml _9_o2 3 
— Kn T^ ^nH')H^ 



50 



^^25? 



H 



-Ml A„2 3 
- 6 + 10^1^2 

= It^'Df^Dyt + \K,lyt + i^gfyt + Ig^yt + ^^ 

I 2 

= \>^DVh>^D + il/fe4 + ^SlVh + ^^2l/b + '^^ 
= In'j^KDyr + |4yr + ^^fl/r + ^^al/r + ^^ 



Au-i 2^ Q^' ^2 9^3 1 33 _4 \ , 11„4\ , 3^2^/ 2 , 3^2^/ 2 , 9gi^o , 3^2„2^/ , 

15g|Kj _ Ml _ ^^2^4 _ X«4_2 _ Jal _n^' 2\ _ 3Kp^ 

4 500 iooy2yi 2oy2yi 4 '^'^o ^^ 2 

/,(!) _ 9^2\ Q„2\ I ^/ ^ _L ^2 I 10X2 I 27gf 9 2„2 , Sfll 

"Ad - -59l^D - ^g2^D + /t£,/tD + /«D + ^^^D + Too" + Io5^2fi'l + T 



k' 2 

2 



1953 ^4 \ I 54^2\2 , 117^2„2\ , p:/l„2\2 51 ^4 \ /i ^/ 2, 

2K'jjKDyr-5K'^Kn\D-^f^'Df^l-Qt^lyt-'^t^lyT-^f^lXD-2K%-78\% + yiK'j^^ + y^K'j^^ + 

H'^'d I 3^2„2^/ I 15g|K'p _ 3537gS _ 1719^2„4 _ 303^4„2 , 291gf _ 3k'p^ _ ^ , 22 _ ^/ 2„,2 
20 "T 2"2yi'^D "T 4 1000 200 "2^/1 40 ^2^1 "•" 8 2 """^-D i/t '^D i/r 

^Id = 6i/b/€D - ffi-^/tD - Qs-Ikd + Qk,d>^h + 2k'^Xd + QK,Dyt + S/t^l/r + S^dAd + 24 + ^ - 

9«2„2 I 9<y| I 9^/ \ 1^/2 
5fl'2fi'l + ^" + "^^D^H + K-D 

,(2) _ 35379? 909^2^4 I 9 ,,2^4 171„,2_4 9,,2_4 , 9K'pgf 1773 „4 , 27 \ „4 , 27 x _4 , 

"•^D - -^00~ + T00fi'2fi'l + To^bS'l ~ ^^tS'l - 2^rfi'l H To~ + -200 '^Dgi + Jo^ogi + iQ^Hgi + 

%9Wi - \^^d9\ - yD9l9l - T9lyl9l - f9bt9f - fgbhl + hWi + f^glt^ogl + 

lylKDgl+^y^tt^Dgl+^yli^Dgl-^glXDgl+^iK'^XDgl+^KD\Dgl-^glXHgl+^K'jj\Hgl+ 

36^ \ ^2 I 291gf 15Kgg| o^/ 3 o^3 1 Q^/ 2„2 9^4„,2 o^/ 2„,2 9^4„,2 o^/ 2„,2 

'^^i^'oyWi - laWr - i^'d^vI + ^gWo - ^vWd - ^ytf^l - 2i/^4 - ^/^d - 44^d - is^d Af, - 
44^1f - 15kd A|^ - ^g^KD - ^yti^D - fytf^D - lytKD-4.K'j^^KD - Qi^'dqIi^d + ^ghli^D + 
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^^9lyl^^D + ^gWtKD + AQglylKD -2lyly^KD + 'ighrf^D + fgl^D - 7k'j,^\d + IS^'^^i^i, - 

18KlXD + 3Qg2KDXD-8K'uKDXD+fgt^H-7K'j^^XH + l8K'j^glXH-l2K'ijy^XH-l2K'oy^XH- 

An'j^ylXu - ISkJjXh + ^t^glnoXn - "i^ylKoXu - 'i^y'^KoXn - Uy'^KoXn - 8k'j^kdXh 
b^l = Qn'^yl + Ak'j^kd + 2k'j,Xd - ^ - 9(?i/s:'^ + ^^glgl + 2«:i,A^ + 2k'^'' + Qn'^^yj + 2/s:'^y2 




A. 5 SM + Majorana Singlet Fermion + Dirac Doublet Fermion 
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A6g 
A6g 



a41 
a4' 



A6il 
A&g 



A6i^ 
A6? 



^ M 



5 



= -i.oro? 



20yiyd 20i'li'M + 50 ^ loyzi/i 



3 



— 4^25^ 4y2y« ^^ 6 ^^ 10^1^2 



9 ~r 9 



2^2 

= i^d'i/'yt + Ig'uVbyt + f^gl&t + Y&gtg?yt + TA~g'uyt + i^i^^'y* + §-^gt 

_9_-^'4 _ 5~/2~/2 _ 9-/23 _ Qg^'^^t _ 9^' 2 3 
lei/d yt 4i/d i/« yi s"d (Jt 16 s"^ ^t 



29 4 I 1 4 

^katyt + 2^2^* 



9 ~r 9 



2 "•" 2 
_ 3 _2~/2„, I 15„2~/2„, I 3_2~/2„, i 15„2~/2„, 1 

- m9i9d yb + iQg29d yb + i^gigu yb + i^g29u yb - 150, 

5 ^'2 2 _ 9-/23 _ '^'g'u^Vb , 5-/22 _ 9-/23 

gi/d ybi/t s9d ilb 16 ^ s9u ybyt s9u (Jb 



;^ty6 + \9tyb - jQg'd^yb - \~g?g'uyb + 



= gd ?^^ I g(.^i/^ 
2 "•" 2 

9 ^' 2„,3 _ 9g'J-yr _ 9 -/ 2„,3 
%9d Ut 16 s"" y-r 

= -g'd' - ^g'd'g'u' + ^g'd'^n - ~g'J + 2g'jXH 

= \grg?^H+'igrg?Xn-^,grg?-Wd'-Mi9rg?+lgr^^^^^ 

WJ - Matg? + lgt>^H + ^glX^ - g - Agl - \g\g\ - %g\gl + ^ + i^ - 
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^t = ^9'uyl - ^ - ^ + ^~9?~9'u + ¥ + ^~9'uy't + ~9'uyl 

9'u ~ ?.9i9uyh + 8 929uyb + '^^9z9uyb 4 yd y«y6 s9u yb ^ 29uybyt 4 ^^ m9i9d 9u + 

39^2-/2-/ I 309 2-/3 , 165^2-/3 , iZ 2~/ 2 , 45 2 -/ 2 , oOrt^rt' 7/2 -L 15„2~/ 2 , 15 2~/ 2 , 
8 929d 9u ^ I609l9u ^ 32 929u ^ 8 9i9uyt ^ 8 929uyt ^ '^^939uyt ^ 8 9l9uyT ^ 8 929uyT ^ 
117gtg(. _ 21g4g; _ 27„2„2~/ _ Qg^i^g; _ 15g;;^g;,3 _ q~/2~/ X _ 21~/2~/„,2 _ 7~/2~/„,2 _ Sg/ _ 

200 4 209l929u 4 4 "-"i/d i/M^H 4 i/d i/ni/t 4yd i/Myr 4 

'^/^'Sx,^ 27~/3„,2 9~/3„,2,M24 27g;g^ 9g(.i/:^ 

41^ = ^g'dv! - ^ - ^ + T + 2^d^«^ + 3^>? + g'dvl 

^ = Wavl + fala'ayl + 20yiy>,^ - 'ic,',hl - 'irnVb + Ig'.vWt - ^ + W^ara? + 
t ^i^:.' + hal-g'^g? + f ^i^:.^;^ + f ^?^>*^ + falg'.y't + 20^1^;,^? + f y?y>^ + f yi^;,!/? + 

ll^gfg^i 21gtg^ 27„2„2~/ ^~9'd^ ^^a'/g'u'^ qn'3X 27~/3„,2 9~/3„,2 ^9'd9'u^ ^n' n' 2 X 

200 4 2oyiy2yd 4 4 "^yrf ^^ 8 yd y* sy^i y^ 4 "^ydyu ^^^ 

21~/~/2„,2 7~/~/2„,2 I 3gd^lf 27g;;i/4 9g;^y4 

4 ydy« y* 4ydy« yr "•" 2 4 4 



A. 6 SM + Majorana Triplet Fermion + Dirac Doublet Fermion 
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= f i/fei/t^d + Wbyt~9l - \yt~glgl + YQglyt~gl - §ytgd - ^y^gl + ^gWgl + ^^?i/t^' 

45,, ^4 _ 27„,3^2 I 165 „2„, ~2 , 29 „4„, , 3„4„, 

i^ytgu -^ytgu + -i^g2ytgu + i^fi'i?/* + 25^2^* 



= lyhgl + Ivbgl 



= fybVtgd - \yb~gl~gl + ^^?i/6^d - %yb~gd - ^yl~gl + W^i?/&^d + ^iybyVgl + ^^?i/6^^ 
iy6^« - ^ylgl + W^2^6^« - ^g\yb + i^2i/b 

= ii/r^i + lyr~gl 

- _3„, ^2-2 , _9_ 2 ~2 _ 45 ~4 _ 27 3~2 , 165^27, n"^ _l _9_rt27, ^^ - 45„, ^4 _ 27 3~2 
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27 ^4-2 I 63„2„2~2 I 4'95 00^2 ~4 153 „4~2 , 9 „2 \ ~2 5 \ ~4 qc \ 2 ~2 , 165 „2 \ ~2 

Too9i9d + To929i9d + ^ ^^^9298 - -^929d + i9i^Hgu - i^ngu - ^^^h9u + -^92^Hgu - 
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|i/fe^d^« + Ivlvlgu - fylgl - fytgu + Iglvlgu + fghlgu + ^Ogiy^gu - XhA + IvtA 
\ylgTgu - f&a - yA + ^gfA + ^gUTgu - ^>^H~gl + l^^Vgu - fyUl - ^yfgu 

17^2 2~ I 45 2 2- ^ 0(\ rP' 1 1"^ 7, — 15,,2^3 _ 9 4~ , 15 2 2~ , 15 2 2~ _ IqI , 87 2~3 

-s9iyt9u + -8 92yt9u + '^^g^ytgu -^yrgu lyrgu + -^giyrgu + Y92yT9u — 2~ ^ 329i9u 

32 92gu ' 200^1^" 12 92gu ' 20^1^2^/" 



4^j = ^y!9d + ^ytgd + ^d^2 + y^ + ^ - ^^?^d - f ^i^d 



C = i^b^/t^^'^ + bbMl - fylgl - ^ytgd + IghVgd + fghlgd + '^^ghlgd - XhMI - ^^^ugl + 
l>^Vgd + lylMl - fytgd - ^yfgd + Y9bt9d + fgbtgd + '^Ogbtgd + lylMl - ^~gTgl - 
\9d~gt + h9l~9d~gl + ^iglMl - ^iyl~gl - \y% + ^ig\yl~gd + ^;glyl~gd - ? + %g\~g\ + 

875 _2~3 I 117^4- 409^4- , 9 ^2_2~ 

"32"5'25'd + 2009l9d - -i2929d + 2o9l929d 

B One-Loop Threshold Corrections 

For a fully self-consistent next-to-leading order calculation, one has to include one-loop threshold 
effects at the weak scale. The pole Higgs and top masses are related to the quartic Xh and top 
Yukawa coupling yt via 

m^(l + 5H) = 2\hv^ (21) 

mt{l + 5t) = ytv, (22) 



where the correction 6h is equal to 

.2 



0£f {xhF, + Fo + xJIF^, + F^p) (23) 



Sh - 

St = St,SM + St,NP (24) 



Here the functions Fi, Fq, F_i are SM contributions to the Higgs quartic which can be found 
in [13], while 5t^sM are SM contributions to the top Yukawa which can be found in [4]. 

The quantities F/vp and Gnp denotes new physics contributions. Corrections for the sin- 
glet/doublet or doublet /triplet fermions may be found in [4], who use a slightly different nota- 
tion. We have calculated the one-loop thresholds for the scalar models utilizing the methods 
of [13, 14]. The corrections to the Higgs quartic from a singlet, doublet and triplet scalars are 
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2 1 2 
= Vflg + /€5f 
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2 1 2 


^?)o = 


= m|, + Kz)f 72 


"^L = 


= m|, + (kd + k'j^)v^/2 



in iH 



Here we have defined ^x = m\/m\. Note that while ms,T,D denote the Lagrangian parameters 
defined in Sec. 2, we have also defined that 

(26) 
(27) 
(28) 
(29) 

which denote the physical masses within the singlet, triplet, and doublet multiplets. 

Finally, we present the threshold correction to the top quark Yukawa coupling from new 
physics, which reads 

al (-^D+ + ^Do + '^xdoXd+ log (^)) sec2 (Ow) 
5tD = -^ ^7 \^^lll (30) 

for the scalar doublet and vanishes in the case of a scalar singlet and triplet. 
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